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What we believe to be a new type of resonant coupling of an incident bulk wave into guided modes 
of a slab with a thick holographic grating is shown to occur in the presence of strong frequency 
detunings of the Bragg condition. This happens through the reflection of the strongly noneigen +1 
diffracted order with the slab-grating boundaries, the resultant reflected waves forming a guided 
slab mode. Rigorous coupled- wave analysis is used for the numerical analysis of the predicted 
resonant effects. Possible applications include enhanced options for the design of multiplexing 
and demultiplexing systems, optical signal-processing devices, optical sensors, and measurement 
techniques. 



I. INTRODUCTION 

Resonances and anomalies of scattering of electromag- 
netic waves in periodic gratings have been extensively an- 
alyzed in the past 0, QitJ ■ This includes resonant gener- 
ation of guided and surface waves by means of bulk waves 
interacting with periodic surface-relief gratings Q, @- 
The associated anomalies of reflection and the strong res- 
onant increase of the scattered wave near the interface 
were called Wood's anomalies 0,0. They occur when 
the scattered wave propagates almost parallel to the in- 
terface of the two media with a grating. The recent rapid 
development of photonics, optical communications, and 
signal processing has resulted in a strong renewed interest 
in these anomalies, mainly due to the possibility of using 
them for the design of multiplexing and demultiplexing 
devices and coupling of bulk electromagnetic waves to 
guided or surface modes 0. All these effects usually 
take place in thin gratings of thickness that is smaller 
than or of the order of the wavelength (e.g., surface-relief 
gratings with small amplitude). 

Recently, several new and unusual resonant effects in 
thick gratings (of thickness that is much larger than 
the wavelength) have been predicted theoretically in 
the geometry of extremely asymmetrical scattering and 
grazing-angle scattering (GAS), i.e., when the scattered 
wave propagates parallel5 or almost parallel to the grat- 
ing boundaries jg. In particular, strong GAS resonances 
have been associated with the generation of a new type of 
mode in gratin gs a nd photonic crystals — grating eigen- 
modes 0T H, EiEl- These modes can be guided by a 
slanted grating alone without conventional guiding effects 



in the structure 0, H, @ ■ At the same time, increasing 
the mean dielectric permittivity inside the grating re- 
gion, i.e., assuming that the grating is formed in a guid- 
ing slab, results in a number of new modes propagating 
in the slab [T^j. These modes have nothing to do with 
the conventional guided slab modes — they are generated 
at wrong angles of scattering and have an unusual field 
structure inside the slab [Io| . 

Grating eigenmodes are formed by all the diffracted or- 
ders interacting with each other inside the grating,and 
thus cannot exist in the absence of the grating U H, 
IToj ]. They appear to be weakly coupled to the bulk in- 
cident wave, which results in strong resonances of scat- 
tering in the grating at angles that do not correspond 
to generation of the conventional guided modes [ll|. In 
fact, grating eigenmodes can exist only if the conven- 
tional guided modes are not effectively generated due to 
large grating amplitude and or grating width, and, vice 
versa, structures in which conventional guided modes are 
effectively generated cannot support grating eigenmodes. 
Therefore grating eigenmodes exist only in sufficiently 
wide gratings with large amplitude (e.g., in photonic 
crystals) . 

Remarkably, grating eigenmodes can exist (and the as- 
sociated resonances can be very strong) in the presence of 
strong frequency detunings of the Bragg condition |9Lll2|. 
In fact, in some structures, the strongest resonances oc- 
cur when the frequency is almost halfway between the 
two Bragg frequencies corresponding to the first- and 
second-order scattering |9|. However, strong resonances 
of scattering at frequencies that are strongly detuned 
from the Bragg frequency normally occur at large grating 
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amplitude (more than w 10% of the mean permittivity 
in the grating) If the grating amplitude is reduced, 
the resonances disappear and eventually are replaced by 
resonances caused by the generation of conventional slab 
modes |To| . 

The aim of this paper is to investigate theoretically 
a new effect that is associated with strong and sharp 
resonances in wide holographic gratings in a guiding slab 
with small grating amplitude and at large (up to « 50%) 
detunings of the Bragg frequency. Remarkably, it will 
be shown that the scattered wave is not only strongly 
noneigen in the slab, but also propagates in a direction 
that does not correspond to a guided mode of the slab. 
Nevertheless, the predicted resonances will be explained 
by generation of the conventional guided modes, but in 
a peculiar way by means of interaction of the noneigen 
(due to strong frequency detunings) scattered wave from 
the slab boundaries. 

II. STRUCTURE AND METHODS 

Consider a slanted holographic grating in a thick di- 
electric slab of width L A (the wavelength in vacuum) 
with the dielectric permittivity given as 

ei for x < 0, 

e(x, y) = e 2 + e g exp(iq • r) + e* exp(iq • r) for < x < L, 
£3 for x > L 

(1) 

where £ 2 is the mean permittivity in the grating re- 
gion; ei and £3 are the permittivities in front of and be- 
hind the grating, respectively; Aei = t\ — e 2 < and 
Ae 2 = £3 — €2 < (guiding slab); e g is the grating am- 
plitude; q = (q x ,q y ) is the reciprocal lattice vector; and 
q = 2n I Lambda, with A as the grating period. The sys- 
tem of coordinates is presented in figure 1; all the media 
are isotropic and nondissipative, and the grating is infi- 
nite in the y and z directions. 

Let a TE electromagnetic wave with amplitude E\q 

1/2 

and wave vector kio (fcio = £1 w/c, where ui is the an- 
gular frequency and c is the speed of light in vacuum) 
be incident onto the grating from the left (x < 0) at an 
angle #10 with respect to the x axis (figure 1). In the 
rigorous coupled-wave theory, the field E(x, y) inside the 
grating is represented in the form of a superposition of an 
infinite number of spatial harmonics (diffracted orders) 
with x-dependent amplitudes E2 n {x) |^|: 
+00 

E(x,y)= ^ E 2n(x) e^q>(ik 2nx x + ik 2ny y - ivt), (2) 

n— — oc 

and the wave vectors k 2 „ are determined by the Floquct 
condition Q: 

k 2 « = k 20 - nq, (3) 

where k 2 o is the wave vector of the incident wave (the 
0th diffracted order) in the grating (slab), with fc 2 o = 

1/2 / 
e 2 ' uj/c. 




< > 

L 

FIG. 1: Thick holographic grating with period A, amplitude 
e g , and reciprocal lattice vector q (q = 2n/A) in a guiding di- 
electric slab of width L. All the materials are assumed to be 
linear, isotropic, and lossless; the dielectric permittivities in 
front and behind the guiding slab are £1 and £3, respectively. 
The mean permittivity in the slab is £2. A bulk TE electro- 
magnetic plane wave with wave vector kio and amplitude Em 
is incident onto the slab at the angle #10 with respect to the 
x axis. The axes of coordinates are also presented. 



Substituting equation (2) into the wave equation in the 
grating gives an infinite set of coupled- wave secondorder 
differential equations that are then truncated and solved 
numerically for the amplitudes of the diffracted orders 
E 2n {x) (see, for example, ,4J). The unknown constants of 
integration are determined from the boundary conditions 
at the grating boundaries. If at some frequency u) = ujq 

1 /2 

and n = 1 we have k 2 \ = fc 2 o = £ 2 <^o/c, then the Bragg 
condition is satisfied for the +1 diffracted order, and u>q is 
called the Bragg frequency (first-order Bragg scattering) . 
In this case, the +1 diffracted order is an eigen bulk wave 
in the material of the slab. If the frequency u> is such 
that k 2 i 7^ k 2 o, then the Bragg condition is detuned by 
the value Aw = lu — luq. 

In this paper we will investigate the frequency response 
of thick holographic gratings (described by equation (1)) 
in a guiding slab. The grating amplitude is assumed to be 
small, so that the GAS resonances do not occur, and res- 
onant generation of conventional slab modes is expected 
in the slab. In Section 3, strong anomalies of scattering 
will be predicted and explained at large (up to w 50%) 
negative frequency detunings of the Bragg condition, i.e., 
when the frequency of the wave uj w luq/2. Rigorous 
coupled-wave theory based on the enhanced T-matrix 
algorithm will be used for the analysis, although in 
some cases the approximate coupled-wave theory based 
on the two- wave approximation |14j will be sufficient. 
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III. NUMERICAL RESULTS 

As an example, consider a structure with the following 
parameters: L — 10 /im; t2 — 5; e g = 2 x 1CP 3 ; #io = 
45°; ei = e 3 = 4.8492, i.e., Aei = Ae 2 = -0.1508 (the 
value of Ae has been chosen so that the critical angle of 
total internal reflection for a wave propagating inside the 
layer is equal to 80°). The Bragg condition is assumed to 
be satisfied precisely at the frequency uq — 1.88 x 10 15 
rad/s (the corresponding Bragg wavelength in vacuum is 
A = l^m). 

Figure 2 presents the dependencies of the relative 
amplitude of the scattered wave (+1 diffracted order) 
\E2i/Eio\ in the grating on the angle of scattering 6>2i- 
This is the angle of propagation of the +1 diffracted or- 
der in the grating with respect to the x axis. The Bragg 
condition is assumed to be satisfied precisely for all con- 
sidered angles of scattering, i.e., the slanting angle and 
the period of the grating are adjusted accordingly for 
each angle of scattering. The dependencies are presented 
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FIG. 2: Dependencies of the relative amplitudes of the +1 
diffracted order I-E21/-E10 on the angle of scattering 621 (i.e., 
angle between the vector lt2i and the x axis in the grating): 
(a) at the front and rear grating boundaries (i.e., at x — 
and x — L) and (b) in the middle of the grating (i.e., at 
x — L/2). The slanting angle for grating fringes and the 
period A are adjusted for each angle of scattering 6*21 so that 

the +1 diffracted order satisfies the Bragg condition precisely 

1/2 

(i.e., tei = e 2 u)a/c). The structural and wave parameters 



are L = 10 (im, Ao (vacuum) = 1 /im, £2 
(the critical angle in the slab is 80°), e g 
angle of incidence 610 = 45° . 



5, ei = e 3 = 4.8492 
= 2 x 1CT 3 , and the 



at the front and rear boundaries of the grating (guiding 
slab) (figure 2(a)) and in the middle of the grating (slab) 
(figure 2(b)). 

Sharp resonances of scattering can be seen in figure 2, 
indicating a strong increase of the scattered wave am- 
plitude at several resonant angles of scattering. At these 
angles, resonant coupling between the incident bulk wave 
and the modes guided by the slab takes place. The max- 
imums in figure 2(b) represent the symmetric modes of 
the slab. This is because the electric field of the anti- 
symmetric modes in the middle of the slab is equal to 
zero (and we do not see the corresponding resonances in 
figure 2(b)). At the same time, at the slab boundaries, 
both symmetric and antisymmetric modes have a nonzero 
electric field. Therefore twice as many maximums can be 
seen in figure 2(a), corresponding to both symmetric and 
antisymmetric modes. 

The patterns of the resonances are almost symmet- 
ric with respect to 621 = 90° in figures 2(a) and 2(b). 
This is because a guided mode in the slab can be repre- 
sented by a bulk wave successively reflecting from the slab 
boundaries |ll| |. Therefore such a mode can equally be 
generated by means of coupling of the incident bulk wave 
into the +1 diffracted order propagating toward the front 
grating boundary (621 > 90°) or toward the rear grating 
boundary (8 2 i < 90°), as long as this diffracted order 
propagates at the angle (with respect to either of the slab 
boundaries) corresponding to a guided mode. Therefore 
the pattern of resonances is symmetric with respect to 
621 = 90°. As a result, the eight maximums in figure 
2(b) correspond to four different symmetric slab modes 
(pairs of maximums that are symmetric with respect to 
6*21 = 90° correspond to the same mode). Similarly, 16 
maximums in figure 2(a) (at the slab boundaries) corre- 
spond to four symmetric and four antisymmetric modes. 
Note, however, that the mentioned symmetry disappears 
when the thickness of the slab becomes of the order of 
or less than the wavelength of the incident wave 0] . In 
this case, the Bragg regime of scattering is replaced by 
the Raman-Nath regime and it is hardly possible to 
speak about a slanting angle of the grating fringes with 
respect to the slab boundaries. This effectively means 
that the reciprocal lattice vector of the grating is parallel 
to the slab, rather than pointing at some angle with re- 
spect to the grating slab boundaries. As a result, the pat- 
tern of the resonances becomes asymmetric with respect 
to 6>2i = 90°, and the maximum(s) at 6> 2 i > 90° p| . 

We emphasize again that thick holographic gratings in 
a guiding slab resonantly couple the incident radiation 
into a slab mode if the angle of propagation of the +1 
diffracted order (for which the Bragg condition is satis- 
fied) is equal to the angle of propagation of the bulk wave 
representing the guide mode of the slab. 

Let us now consider the frequency response of scatter- 
ing in the thick holographic gratings in a guiding slab. 
For example, consider the same structure with the thick 
holographic structure in a slab of thickness L = 10 ^m, 
with e 2 = 5; e g = 2 x 10~ 3 ; 6 10 = 45°; ei = e 3 = 4.8492, 
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i.e., Aei = Ae2 = Ae = —0.1508. However, this time 
we fix the orientation and the period of the grating so 
that the Bragg condition for the +1 diffracted order is 
satisfied precisely at lo — loq and the angle of scattering 
is 02i = 70°. At the mentioned structural parameters, 
we change the frequency of the incident wave. The de- 
pendencies of the amplitude of the +1 diffracted order on 
frequency detuning of the Bragg condition Auj are pre- 
sented in figures 3(a) and 3(b) at the boundaries and in 
the middle of the slab, respectively. 
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FIG. 3: Dependencies of the relative amplitudes of the +1 
diffracted order |i?2i/-Eio| on frequency detuning Alo: (a) at 
the front and rear grating boundaries (i.e., at x = Oand x — L) 
and (b) in the middle of the grating (i.e., at x = L/2). The 
structural parameters are the same as in figure 2. However, 
this time we fix the grating orientation and period so that the 
+ 1 diffracted order satisfies the Bragg condition at #21 = 70° 
and u! — luo = 1.88 x 10 15 rad/s (Ao = 1 /im). 

The main unusual features of these dependencies (fig- 
ure 3) are the strong and sharp resonances of scatter- 
ing at the frequency detunings Aw that are just w 5-6 
times smaller than the Bragg frequency uio — 1.88 x 10 15 
rad/s. Note that it is hardly possible to expect any 
scattering in the grating at the considered small grat- 
ing amplitude e g and such large frequency detunings. It 
is obvious that these resonances must be related to the 
resonant generation of some type of localized or guided 
modes in the structure. Further decrease of the grat- 
ing amplitude results in a rapid increase of the height 
and sharpness of the considered resonances. However, 
the angles at which these resonances occur are hardly 
affected by reducing the grating amplitude. If the grat- 
ing amplitude tends to zero, then the resonances tend to 
infinity, which means that the corresponding modes still 



exist, but they become uncoupled to the incident wave. 
This demonstrates that these modes can exist in the ab- 
sence of the grating. Therefore they cannot be a type of 
grating eigenmode 0, H, El Ell but must simply be the 
conventional guided modes of the dielectric slab. This 
has also been confirmed by considering the field distribu- 
tions in the grating and slab at the resonant frequency 
detunings. It can be shown that these distributions are 
practically identical to those corresponding to the con- 
ventional guided modes in the slab. As a result, it can 
be seen that the resonances in figure 3(b) correspond to 
the generation of the symmetric modes of the slab, and 
the additional maximums in figure 3(a) correspond to the 
antisymmetric modes. In particular, the leftmost peaks 
in figures 3(a) and 3(b) (at the largest Acu) correspond to 
the fundamental guided mode of the slab, and the peaks 
at smaller detunings represent the higher-order guided 
modes. 

However, there is still a question as to how these slab 
modes could be generated by the incident wave at the 
mentioned strong (up to « 20%) frequency detunings 
Ao; of the Bragg condition, when the +1 diffracted order 
in the grating is expected to be a strongly noneigenwave 

1 /2 

with the wrong wave vector k2i = k2o q, fei 7^ e 2 w / c 
(u> = ujq + Alu), which does not correspond to a bulk wave 
propagating in the material of the slab. Moreover, it can 
be shown that at each of the resonant detunings (figure 3) 
the angle between the vector k 2 i and the slab boundaries 
is significantly different from the angle that one would 
expect for a guided slab mode. Therefore noneigen oscil- 
lations of the electromagnetic field corresponding to the 
+1 order should be driven by the incident wave in the 
grating and must rapidly tend to zero if e g — > 0. Such 
a noneigenwave in the slab cannot form a guided mode. 
This, however, is in obvious contradiction with the ob- 
tained rigorous numerical results (figure 3). 

To understand the physical origins of the resonances 
observed in figure 3, consider figure 4. If Auj = (i.e., 
lu = loq), then the wave vector of the incident wave inside 
the grating and slab is k2oo- In this case, we assume 
that the Bragg condition is satisfied and thus the +1 
diffracted order is a propagating eigenwave in the slab. 
Thus the wave vector of this wave k2io = k2oo — Q is equal 

1 /2 

in magnitude to k2oo, i.e., ^210 = ^2 w o/ c (the angle 
between the k2io vector and the x axis is #210 = 70° — 
see above). 

Introducing nonzero detunings of the Bragg condition 
by reducing the frequency of the incident wave by Aw 
results in reduction of the wave vector of the incident 
wave k2o (figure 4). At all frequency detunings, the wave 
vector k 2 i of the +1 diffracted order is given by the Flo- 
quet condition (equation (3)) with n = 1. Therefore it 
can also be seen that if Auj < 0, then | k2 1 1 > |k2o|, and 
the angle between the k2i vector and the x axis increases 
with increasing magnitude of the detuning (figure 4). 

As mentioned above, the wave vector k2i has a wrong 

1/2 

magnitude (kzi 7^ e 2 w/c) and direction (that does 
not correspond to a guided mode in the slab). The 
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FIG. 4: Scheme of coupling between the bulk waves and 
slab modes by means of the interaction of the noneigen +1 
diffracted order with the grating boundaries: k2oo and k2io 
are the wave vectors of the Oth and +1 orders in the slab when 
the Bragg condition is satisfied (Aw = and w = wo); k2o and 
k2i are the wave vectors of the Oth and +1 orders at nonzero 
frequency detuning Aw < 0; k2i r and k 2 i r are the wave vec- 
tors of the waves resulting from the reflection of the noneigen 
+1 diffracted order (with the amplitude E21S and the wave 
vector k2i) from the grating boundaries: k2i ry = fc 2 i ry = k2i y 
(Snell's law) and k^ir = k' 21r = (wo + Aw)(e2) 1//2 / c (i.e., the 
reflected waves are eigen for the slab). When the angle of 
propagation of these reflected waves 6 r corresponds to one of 
the guided slab modes, resonant generation of this mode oc- 
curs with a strong resonant increase of the amplitude of the 
+ 1 diffracted order in the grating (figures 3, 5, and 6). This 
is because the reflected waves are formally included in the +1 
diffracted order (see equations (2) and (5)). 



noneigenwave with this wave vector interacts with the 
slab-grating boundaries and generates a periodic field 
pattern along these boundaries. According to the Huy- 
gens principle, this pattern will generate a reflected bulk 
wave in the slab with the wave vector k 2 i r (figure 4). 
The tangential (to the boundaries) component of this 
vector k 2 i r is equal to the y component of the wave vec- 
tor k 2 i: ki\ry — k 2 i y (Snell's law). However, since this 
reflected wave is not driven by scattering (its wave vector 
does not satisfy the Floquet condition of equation (3)), 

it must be either an evanescent wave in the grating (if 
1/2 

k2iry > e 2 w / c ) or a propagating eigenwave in the ma- 

1/2 

terial of the slab (if k2i ry < e 2 LJ / C )- Therefore the x 
component of its wave vector k 2 i r is given as 



k 2 irx = -Je 2 uJ 2 /c 2 - k\ Xy ^ -k 2 i x . (4) 



If this reflected wave is an eigenwave in the slab (i.e., 

1 /2 

k.2\ry < e 2 w /c = &2ir) an d it propagates at an angle 
that corresponds to one of the guided modes in the slab, 
then this guided mode will be resonantly generated in 
the slab (figure 3). Thus the guided modes of the slab 
in this case are generated by means of the interaction of 
the noneigen +1 diffracted order with the slab bound- 
aries. Only the presence of the boundaries (structural 
discontinuities along the x axis) results in transforma- 
tion of the wrong wave vector of the +1 diffracted order 
in the presence of strong frequency detunings into the 
correct vector corresponding to a guided slab mode. 

Since changing the frequency results in changing the 
direction of the k 2 i vector and its magnitude, the cor- 
responding direction of the wave vector of the reflected 
eigenwave in the slab k 2 i r will also change. Thus, ad- 
justing frequency detuning, we automatically adjust the 
direction and magnitude of the vector k 2 i r . When the 
detuning is such that k 2 i r corresponds to a guided slab 
mode, resonant generation of such a mode occurs (fig- 
ure 3). The resonant detuning must be strong so as to 
result in a significant adjustment of the magnitude and 
direction of the k 2 i r vector. Since the detuning is strong, 
the amplitude of the noneigen +1 diffracted order in the 
grating is very small. This means that the coupling of 
the incident wave into the corresponding slab mode is 
very weak, which results in very strong and sharp res- 
onances in figure 3. Decreasing the grating amplitude 
results in a further decrease of the coupling and the cor- 
responding increase of the height and sharpness of the 
resonances (figure 3). However, the resonant frequency 
detunings will hardly depend on the grating amplitude, 
since it does not affect the directions and the magnitudes 
of the considered wave vectors (figure 4). This is clearly 
confirmed by the numerical results. 

The height of the observed resonances (figure 3) 
is practically infeasible — it will be hardly possible to 
achieve such resonances experimentally because this 
would require very large relaxation times (due to weak 
coupling). This example has been used to highlight the 
considered new type of coupling and to demonstrate that 
the resulting resonances are associated with the conven- 
tional guided modes generated in an unusual way by 
means of reflection of the noneigen +1 diffracted or- 
der from the slab boundaries. A more reasonable sit- 
uation occurs when the grating amplitude is increased 
(e g = 0.08) and the corresponding resonances are sig- 
nificantly lower, as shown in figure 5. As can be seen 
in figure 5, the positions of the resonances are practi- 
cally the same as in figure 3 with only w 1% reduction 
of the resonant frequency detunings. This demonstrates 
very weak dependence of the resonant detunings on grat- 
ing amplitude. This is expected, since only large grating 
amplitudes are expected to noticeably change the guiding 
properties of the slab and the corresponding wave vectors 
of the slab modes. 

The fact that the resonances in figures 5(a) and 5(b) 
still correspond to conventional guided modes is again 
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FIG. 5: Dependencies of the relative amplitudes of the +1 
diffracted order -E21/-E10I on frequency detuning Alu: (a) at 
the front and rear grating boundaries (i.e., at x = and 
x — L) and (b) in the middle of the grating (i.e., at x = L/2). 
The structural and wave parameters are the same as in figure 
3, except for the increased grating amplitude: e 9 = 0.08. 
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FIG. 6: Dependencies on x of the relative amplitudes of the 
+1 diffracted order I-E21/-E10I inside the grating and slab for 
the three leftmost resonances in figure 5(a): (a) e g = 0.08, 
(b) e g — 0.1. Curves 1, 2, and 3 in (a) correspond to the 
leftmost, second leftmost, and third leftmost resonances in 
figure 5(a). Resonant detunings are Awi ~ —3.71 x 10 14 
rad/s, Auj 2 » -3.67 x 10 14 rad/s, Alo 3 « -3.59 x 10 14 rad/s 
for curves 1-3, respectively, in both (a) and (b). 



confirmed by the analysis of the field distribution in the 
+1 diffracted order inside the slab at each of the res- 
onant frequency detunings. For example, the obtained 
field distributions inside that slab for the three leftmost 
maximums in figure 5(a) are presented in figure 6(a). 
It can be seen that these distributions are very close to 
those in the first three conventional modes guided by the 
slab. Therefore the leftmost maximum in figure 5(a) cor- 
responds to curve 1 in figure 6(a) (fundamental mode 
of the slab), the second- from-the- left maximum in fig- 
ure 5(a) corresponds to curve 2 in figure 6(a) (the first 
antisymmetric mode), etc. 

At the same time, it can be noticed that curve 1 in 
figure 6(a) displays small bumps, which is not a typical 
feature of the field distribution corresponding to the fun- 
damental slab mode. These bumps become more pro- 
nounced with increasing grating amplitude (see figure 
6(b) for € g — 0.1). Curve 2 in figure 6(b) also displays 
similar irregularities (unlike curve 2 in figure 6(a) for 
e g = 0.08). This is related to the effect of the grating 
on the guided modes. This effect rapidly increases with 
increasing grating amplitude (compare figures 6(a) and 
6(b)). Eventually, when the grating amplitude becomes 
sufficiently large, the structures of the guided modes sub- 
stantially change,and they may be transformed into grat- 
ing eigenmodes 0, 0, 0> El- Note also that the effect 
of the grating is more pronounced for lower slab modes 



(compare the curves in figures 6(a) and 6(b)). 

It is also important to note that the reflected wave 
with the wave vector k2i,- does not appear explicitly in 
the coupled- wave expansion of equation (2). Rather, it is 
formally included in the amplitude of the +1 diffracted 
order. Indeed, the +1 diffracted order in equation (2) 
can be represented as 

E(21)(x) exp(ik 2 i • r — iut) = 
[E(21s)(x) + E ( 21r) exp(ik 2 i rx - ik 2 i x ) 

+E'(21r) exp(-ifc 2 irz - ifeirc) 

x exp(ik2i • r — 'nut), (5) 

where E2\ s {x) is the x-dependent amplitude of the +1 
diffracted order resulting from the direct interaction 
(scattering) of the incident wave with the grating, and 
the other two terms in the square brackets are due to the 
two waves reflected from the rear grating boundary [Ei\ r 
with k,2\rx < 0; see equation (4)) and the front grating 
boundary (E 2 i r ). All three waves have the same y com- 
ponents of their wave vectors (Snell's law, or conserva- 
tion of the tangential component of momentum) , whereas 
the x components are significantly different: k 2 i x , A^irm 
and — k2irx, respectively. Because the y components of 
the wave vectors are the same for all three waves, these 
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waves can formally be included in the + 1 diffracted order 
in the expansion in equation (2), which eventually results 
in a strong resonant increase of the overall amplitude of 
this order at the resonant frequency dctunings (figure 3). 
Actually, these are the amplitudes E%x r and E' 2 \ r that 
experience a strong resonant increase at the resonant de- 
tunings, whereas the amplitude Ei\ a {x) remains small (it 
tends to zero as e g — > 0). In this case, because the am- 
plitudes of the reflected waves become predominant at 
the resonant detunings, the field distribution inside the 
grating and slab indeed appears to be almost identical to 
that in the corresponding slab modes (figure 6). 

It is possible to derive analytical equations for the reso- 
nant detunings in a slab with a thick holographic grating. 
For this purpose, we will need to determine the magni- 
tude and components of the wave vectors k2i and k2i r . 
From figure 4, it can be seen that the magnitude of the 
wave vector k2i is determined by the equation 

^21 = ^210 + ^^20 _ 2fc 2 ioAfc 2 o cos(2a), (6) 

where a is the angle between k 2 oo and the grating fringes 
(figure 4) and 

Auj 

Afc 20 = fc 200 - fc 20 = - — V^i- (?) 

It follows that, since 2a is the angle between k2oo and 
k 2 i (figure 4), 

2a — 9 2 io — 6*200- (8) 

It can also be seen that the angle ft between the two 
vectors k2i and k2o is determined by the equation 

q 2 = A|i + &20 - 2fc 2 ifc 20 cos(/3). (9) 

Then the angle # 2 1 between the vector k2i and the x axis 
is 

021 = ^200 + P. (10) 

Using equations (6)-(10), we find the y component of the 
wave vector k2i, and thus the components of the k 2 i r 
vector: 

foly = hi Sin 6>21 = k 2 lry, (11) 
k 2 lrx = \Jt2UJ 2 /c 2 - fcf-L Sin 2 6*21 (12) 

From here, the angle of propagation of the reflected wave 
9 r (figure 4) is 

9 r = arctan ( ^) . (13) 

For example, consider the leftmost peak in figure 
3. It corresponds to the frequency detuning Au>i w 
-3.7279 x 10 14 rad/s. Using equations (6)-(13), we 
obtain 6 r w 88.55°. The conventional analysis of the 



guided modes [llj in the considered slab at the frequency 
lu = uj q + Auj 1 demonstrates that this is exactly the angle 
of the wave propagation in the slab, which corresponds to 
the fundamental slab mode. Thus the leftmost maximum 
in figure 3 corresponds to the resonant generation of the 
fundamental slab mode at the frequency lu — ujq + Au>i 
in an unusual way through the noneigen +1 diffracted 
order, interacting with the slab boundaries. Using a sim- 
ilar procedure, it can be shown that the other resonances 
in figure 3 are associated with the generation of higher 
slab modes by means of the same mechanism. 

It is important to note that, as can be seen from figure 
4, an increase in the magnitude of the frequency detun- 
ing Aw (i.e., a reduction of fc 2 o) results in an increase 
in the angle 6 r . This is the reason why the resonance 
associated with the generation of the lower-order guided 
modes (corresponding to larger angles 6 r ) is found at 
larger detunings (figures 3 and 5). 

It is also important to emphasize that the considered 
effect of resonant coupling at strong resonant frequency 
detuning occurs in a wide range of angles of scattering 
and incidence. The angle 6*210 = 70° has been chosen 
just as an example, and a similar pattern occurs for other 
scattering angles. However, the actual values of the res- 
onant frequency detunings will depend on the angles of 
scattering and incidence. 

Increasing the grating-slab width and/or the mean 
permittivity in the slab results in an increase in the num- 
ber of guided modes. Therefore this will also result in an 
increase in the number of different resonant frequency 
detunings. 

It can be seen that figure 3 displays only seven reso- 
nances, whereas in figure 2 we have eight different maxi- 
mums corresponding to eight different modes supported 
by the slab. It can also be shown that the number of 
resonances at strong frequency detunings (Aw < 0) de- 
creases with decreasing angle of scattering B\ and increas- 
ing I A10 1 . This is because such strong negative frequency 
detunings result in significant reduction of the frequency 
of the guided slab modes. Therefore the thickness of the 
slab is effectively reduced compared to the wavelength, 
and the number of modes that can be sustained by the 
slab decreases. For example, it can be seen that this is 
the highest antisymmetric mode (shown in figure 2) that 
does not result in a maximum (i.e., does not exist) in 
figures 3 and 5. 

It can be seen that the considered new way of cou- 
pling between an incident electromagnetic wave and slab 
modes in the presence of a thick holographic grating can 
occur not only by means of the +1 diffracted order but 
also by means of any other diffraction order in the grat- 
ing (e.g., +2, —1 orders, etc.). In this case, for example, 
the strongly noneigen +2 diffraction order in the grat- 

1 /2 

ing (fc 22 7^ e 2 w / c ) interacts with the slab boundaries 
and, if the detuning and the direction of the vector k22 
are appropriate, the resultant reflected wave may form a 
guided mode. An equation similar to equation (13) can 
also be derived in this case. However, because the higher 
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diffracted orders are much weaker than the +1 order in 
gratings with small amplitudes e g , the efficiency of this 
coupling will be much lower, and the corresponding res- 
onances will be very sharp. This makes these resonances 
less achievable for the real structures, because they will 
correspond to large relaxation times. This problem could 
be overcome by increasing grating amplitude, but this 
will result in a noticeable effect of the grating on the 
resultant guided modes. Therefore the described new 
type of coupling of bulk waves and guided slab modes 
can practically be achieved mainly for the +1 diffracted 
order (figures 5 and 6), although the principles of this 
approach can be applicable to any other higher order of 
scattering. 

IV. CONCLUSIONS 

In this paper we have described a new method of cou- 
pling an incident electromagnetic wave into a guided 
mode of a thick slab with a holographic grating. As a 
result, strong resonances were shown to exist at very 
large negative frequency detunings of the Bragg condi- 
tion in the grating. These resonances rapidly increase 
with decreasing grating amplitude, despite the strong 
(up to w 20%) frequency detuning of the Bragg condi- 
tion. It was demonstrated that it is the reflection of the 
strongly noneigcn +1 diffracted order from the slab grat- 
ing boundaries that may result in a reflected eigenwave 
that forms a conventional guided mode of the slab. Thus 
a new type of resonant generation of the conventional 
slab modes in a thick slab with a holographic grating 



was demonstrated to occur by means of the intermediate 
strongly noneigen +1 diffracted order interacting with 
the slab-grating boundary. 

The predicted resonant coupling of the incident wave 
into a thick slab waveguide with a holographic grating 
may also occur by means of an nth diffracted order (e.g., 
n = +2, n = —1, etc.). In this case, if the grating ampli- 
tude is small, the coupling efficiency should significantly 
reduce with an increase in the intermediate diffracted or- 
der because the amplitude of this order rapidly decreases 
with increasing n. This results in a rapid increase of the 
height and sharpness of the corresponding resonance with 
increasing n. 

Rigorous coupled-wave theory was used for the analy- 
sis of scattering, and the example of bulk TE incident 
electromagnetic waves was considered. However, it is 
important to note that if the grating amplitude is suffi- 
ciently small, the approximate coupled-wave theory Q 
may also be used for the approximate determination of 
the predicted resonances at n = 1, if reflections at the 
slab boundaries are taken into account (see, for exam- 
ple, ^|). At the same time, for other values of n, the 
rigorous coupled- wave theory is essential even if the grat- 
ing amplitude is small. This is because the approximate 
theory based on the two-wave approximation [14j simply 
excludes the higher orders that are essential intermedi- 
aries of the resonant coupling. 

Possible applications of the discovered new type of cou- 
pling include enhanced options for multiplexing and de- 
multiplexing, resonant excitation of slab waveguides, de- 
sign of new optical sensors, and measurement techniques. 
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